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Abstract 

A nonthermal particle acceleration mechanism involving the interaction of a charged particle 
with multiple magnetic islands is proposed. The original Fermi acceleration model, which assumes 
randomly distributed magnetic clouds moving at random velocity V c in the interstellar medium, 
is known to be of second-order acceleration of 0(V c /c) 2 owing to the combination of head-on 
and head-tail collisions. In this letter, we reconsider the original Fermi model by introducing 
multiple magnetic islands during reconnection instead of magnetic clouds. We discuss that the 
energetic particles have a tendency to be distributed outside the magnetic islands, and they mainly 
interact with reconnection outflow jets. As a result, the acceleration efficiency becomes first-order 
of 0(Va/c), where Va and c are the Alfven velocity and the speed of light, respectively. 
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Particle acceleration during magnetic reconnection is a fundamental problem in colli- 
sionless plasma, and it is now widely recognized that magnetic reconnection is one of the 
important nonthermal particle acceleration processes in our universe. The nonthermal parti- 
cles observed in solar flares and substorms in the Earth's magnetotail are known to be driven 
by reconnection , and magnetic reconnection now receives a great deal of interest in 
many astrophysical settings such as pulsars, magnetars, galaxy clusters, and active galactic 
nucleus jets 3|, |4|. 

There has been much acclaimed research into particle acceleration through reconnection. 
The most elementary process is meandering/Speiser acceleration in the diffusion region with 
a weak magnetic field js, 6]. The charged particles inside the diffusion region can be directly 
accelerated almost along the reconnection electric field until they escape the diffusion region. 
However, the diffusion region is very small compared with the total system of reconnection, 
and the amount of energetic particle flux produced by this process is thus not necessarily 
enough to explain the observation of nonthermal energy density, which seems to account for 
at least several tens of percent of the total energy density 

Another mechanism is betatron-type acceleration around the pileup magnetic field re- 
gion where Alfvenic outflow jets in the reconnecting magnetic field collide with preexisting 
plasma at rest 7|, |8j. In the pileup region, the plasma can gain energy proportional to the 
magnitude of the compressed magnetic field. In addition to the betatron-type acceleration, 
if the gyroradius of a charged particle is of the order of the magnetic field curvature ra- 
dius, the chaotic particle behavior can enhance the particle energization Additionally, 
the electromagnetic waves generated by reconnection can be an agent of the nonadiabatic 
process [7j. 

Another important acceleration process can occur for reconnection with multiple magnetic 
islands. For a long current sheet, many elongated magnetic islands can form from tearing 
or plasmoid instability [10]. An elongated magnetic island then starts to shrink, and the 
particles trapped inside the island can be accelerated through a Fermi- type process If 
the gyroradius of the accelerated particles exceeds the size of the island, the particles escape 
the island and their energization ceases. 

Turbulent magnetic reconnection has been also discussed as a possible mechanism of 



particle acceleration 



121 ] . The magnetic tension force of a tangled magnetic field line in 



turbulent reconnection is suggested to be as an important agent of particle acceleration [13 ]. 
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FIG. 1: Time evolution of multiple magnetic islands. Panels (a-c) show the plasma density (color 
contour) and magnetic field lines (white), and panel d shows the energy spectrum at t/r c = 320. 
The dashed line is the Maxwellian fit. Panels (e and f) show the densities for the mid-energy and 
high-energy particles, respectively. 

However, the acceleration efficiency of this process remains an open question. 

In this letter, we investigate particle energization in the interaction of a charge particle 
with many magnetic islands, and propose a new acceleration process extended from the 
original Fermi acceleration mechanism [jjj]. In the original Fermi acceleration mechanism, 
particles gain energy stochastically during head-on and head-tail collisions of particles with 
magnetic clouds or mirrors as the scattering objects. During the stochastic acceleration, the 
increase in particle energy is known to be second order of V c /c (V c and c are the velocity 
of the random motion of the magnetic cloud and the speed of light, respectively), and 
is apparently a slow process. Here we study the dynamic evolution of many magnetic 
reconnection sites instead of magnetic clouds, and we show that the energetic particles have 
a tendency to interact with the reconnection outflow jets. The acceleration efficiency is thus 
strongly enhanced relative to that of the original Fermi acceleration, and the energy increase 
becomes first order of Va/c, where Va is the Alfven velocity. 

To study the stochastic acceleration during the interaction with magnetic islands, simula- 
tions are carried out with a relativistic particle-in-cell (PIC) code. Shown in Figure 1 is the 
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time evolution of the magnetic islands and the energy spectrum. As the initial condition, 
four Harris current sheets of the relativistic pair plasmas with peak density uq superim- 
posed on uniform background density n& = 0.05n are set up in two-dimensional space of 
xxy = 512x512 girds. The initial plasma temperature is assumed to be T / mc 2 = 5 for both 
electron and positron, and the drift velocities of electrons and positrons are v^/c = ±0.5, 
respectively. The plasma skin depth c/Q p = 1.12 A, the gyroradius c/Q c = 0.25A, and 
the magnetization parameter (Q c /Q p ) 2 = 20, where the grid size A. The total number of 
particles is 512 x 512 x 200. The double-periodic boundary condition is assumed for both 
x and y directions. A set of Maxwell's equations are solved in the Fourier space so that 
the numerical dispersion error can be suppressed. The simulation results are presented in 
normalized units: the magnetic field is normalized to the asymptotic value of the antiparal- 
lel component of the Harris model Bo, the density to no, length to the initial plasma sheet 
thickness A = 5A, and time to the light transit time r c = A/c. A uniform guide magnetic 
field B z = Bq is added. 

Shown in Figure la-c is the time evolution of the plasma density (color contour) and 
the magnetic field lines (white). In the early stage of the linear tearing mode instability 
{t/ T c — 80 in panel a), approximately eight magnetic islands form in each current sheet, the 
time evolution of which is consistent with our previous simulation As time goes on (b 
and c), the magnetic islands grow by coalescing and merging, and the islands from adjacent 
current layers overlap. At t/r c = 320 (panel c), two large magnetic islands containing most 
of the plasma density are seen. 

During the above evolution, many nonthermal particles are generated. Figure Id shows 
the energy spectrum, where the horizontal and vertical axes are the Lorentz factor 7 of the 
particle energy and the number density Ni^y). The solid line is the energy spectrum at 
t/r c = 320, and the dashed line is a Maxwellian spectrum fitted by T/mc 2 = 7.3. During 
the evolution of reconnection, the initial thermal plasmas are heated and accelerated. There 
is clear formation of nonthermal particles above 7 > 40, and the maximum energy can reach 
7 ~ 10 3 , for which the gyroradius is almost half the size of the simulation system. The 
nonthermal particle occupies about 66 % of the total particle energy. 

How and where energetic particles are accelerated is now investigated. We plot the 
densities for the energetic particles in panels e and f of Figure 1. We calculate the densi- 
ties by integrating the velocity distribution function with the particle momentum accord- 
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FIG. 2: Illustration of (a) the plasma sheet reconnection in the early stage, (b) the reconnec- 
tion during magnetic island merging/coalescence in the late stage, and (c) the interaction of an 
energetic particle with magnetic islands. Black lines show magnetic field lines, red arrows show 
the reconnection outflows with Alfven speed Va, and the blue dashed line is the trajectory of the 
energetic particle. 

ing to f Pmax f(p)dp. Panels e and f show the densities of the energetic populations with 

" Pmin 

{Pmin,Pmax) = (50,200) and (200,oo), respectively. {{p m in,Pmax) = (0, oo) for panels a-c. 
) The gyroradius for the energetic particle with p/mc = 200 is about 10A, which roughly 
corresponds to the size of the magnetic island during the island merging stage. In contrast 
to the total plasma density, we see that (1) the middle energy range (panel e) is localized 
around the magnetic island, and (2) the high-energy plasma (panel f) forms a void structure 
inside the islands. We find that the thermal plasma is mainly confined in the islands (i.e., 
the weak magnetic field region), while the energetic particles are preferentially distributed 
outside the islands (i.e., the stronger magnetic field region). 

The reason why energetic particles have the tendency to be localized in the strong mag- 
netic field is simply illustrated in Figure 2. Reddish regions surrounded by black lines are 
the magnetic islands and red arrows are the Alfvenic reconnection jets. In the early evolu- 
tion of magnetic reconnection (panel a), the energetic particles are known to be generated 
in and around the X-type reconnection region, and they are ejected into the plasma sheet 
in association with the reconnection outflow jet jf]]. On the other hand, in the later phase 
of magnetic coalescence/island merging (panel bj, the reconnection outflows are always di- 
rectly towards the strong magnetic field region [15J. That is, the thermal plasma is supplied 
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FIG. 3: A typical particle trajectory during multiple magnetic reconnections. The particle ener- 
gization history (a), and the particle trajectories superposed on the magnetic field structures at 
tjr c = 180 (b) and at t/r c = 200 (c). Color contours show the magnitude of the magnetic field 
(B 2 + By) 1 / 2 , and black arrows and thin white lines are the plasma flow vector and the magnetic 
field lines in the reconnection plane. 

by the two merging magnetic islands, and energetic particles are generated in the current 
sheet sandwiched by two magnetic islands and are ejected into the strong magnetic field 
region. The nature of particle energization during the merging of the islands explains the 
localization of energetic parties in the strong magnetic field region seen in Figure If. 

To study the mechanism of particle acceleration in the strong magnetic field region in 
detail, we plot a typical particle trajectory and its energy history as a function of the 
magnetic intensity in Figure 3. Figure 3a shows the energy history from the initial state at 
t s = to the end of the simulation run at t e = 320r c . The horizontal and vertical axes are 
the magnitude of the magnetic field (B xy = (B 2 + By) 1 / 2 ) and the particle energy (7 — 1) 
normalized by the rest-mass energy, respectively. At the initial time t = 0, the particle is 
located around (B xy , e) = (1, 30), which is denoted t = t s . As time passes, the particle gains 
energy by moving around the reconnection region. Note that, roughly speaking, B xy < 0.4 
is inside the magnetic island, while B xy > 0.6 is outside the magnetic island and within 
the strong magnetic field. In the early evolution, the particle is accelerated in the relatively 
weak magnetic field regions with B xy < 0.3, and we checked that the acceleration occurred 
in and around the magnetic diffusion region. This is meandering/Speiser acceleration js]. 
After the initial acceleration in and around the diffusion region, the accelerated particle is 
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ejected from the plasma sheet into the strong magnetic field region during the island merging 
process, and it moves around the strong magnetic field. After a while, an accelerated particle 
encounters another magnetic reconnection region in the island merging stage. 

Figures 3b and 3c show the interaction of a charged energetic particle with the reconnec- 
tion jet of randomly moving magnetic islands. The center panel is the particle trajectory 
superposed on the magnetic field line at to = 180r c , and the right-hand panel is that at 
t\ = 200r c . We see that the charged particle can gain energy during the interaction of the 
charged particle and the reconnection jet with a relatively strong magnetic field of B xy ~ 0.8. 
After t > t\ the particle can further gain energy in the strong magnetic field region. The tra- 
jectory shows clearly that particle energization occurs during the interaction of the particle 
and the magnetized plasma flow of the reconnection region. 

Based on the above simulation result, we propose a new particle acceleration mechanism 
as shown in Figure 2c. In contrast with the original Fermi acceleration mechanism as a 
stochastic means by which charged particles collide with magnetic clouds or mirrors [jjj], 
we introduce magnetic islands instead of magnetic clouds or mirrors. During the merging 
of magnetic islands, energetic particles have the tendency to be distributed outside the 
magnetic islands, and the particles collide stochastically with magnetic islands from the 
outside. Those particles interacting with the reconnection outflow jets with Alfven velocity 
Va will gain energy, and those particles involved in head-tail collisions with the magnetic 
islands may lose energy. However, the random velocity V c may be lower than the thermal 
velocity V s . If V c > V s , the plasma will be heated by a shock wave until V c < V s is satisfied. 
Therefore, for an active reconnection region where the magnetic field energy dominates the 
system (i.e., the plasma (3 < 1), we may assume V c < V s < Va, and the energy loss due 
to the head-tail collision may be negligible. The detailed physics of the random velocity 
remains an open issue. 

In our simple particle acceleration model, we did not specify the scattering process in 
detail, because if we could find a Lorentz transformation into the de Hoffman-Taylor frame 
during the scattering of a particle, the particle energy gain does not depend on the local 
scattering process. If the particles can be reflected back in Alfvenic jets, the particle energy 
gain (Se) is proportional to the Alfven speed; i.e., 8s /s oc Va/c. 

In this letter, we presented only the result for the relativistic plasma case, but we con- 
firmed that nonrelativistic reconnection has also the same behavior. The reconnection out- 
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flow velocity does not depend on the relativistic or nonrelativistic plasmas, and the outflow 
always has the Alfven speed Va- Once energetic particles are ejected from magnetic islands, 
our stochastic reconnection process performs well. The efficiency of energetic particle ejec- 
tion from the islands (i.e., the so-called injection problem) is an important consideration of 
this process, and it seems to depend on the plasma dynamics inside the merging current 
sheet. 

The spatial size of the simulations presented here is small compared with any natural sys- 
tem, but we also confirmed in another simulation run that maximum attainable energy was 
about four times larger for a four times larger simulation box. The maximum attainable 
energy e max can be simply expressed by ■y max = £max/ m c 2 = eBL/2mc 2 ~ 10 3 (L/100A), 
where L is the simulation box size. This estimation agrees well with Figure If. This mech- 
anism naturally predicts a power law energy spectrum by taking account of particle loss 
from the system, but the power law index may depend on the dynamics of the multiple 
reconnection and the process of particle loss from the active region of island merging. To 
save the computational time, pair plasma has been assumed, but it would be interesting to 
compare the particle acceleration efficiency for the proton and electron in future study. 

Particle acceleration in multiple magnetic islands is important in many astrophysical 
settings. For example, the pulsar wind is known to form a striped wind structure of the 
toroidal magnetic field in the equatorial plane, and it has been proposed that magnetic 
reconnection is an important acceleration process that explains synchrotron radiation from 
the Crab nebula 4J • At the boundary between the pulsar wind and the synchrotron nebula, a 
relativistic perpendicular shock is believed to form [16j . In addition to the shock acceleration, 
the generation of multiple magnetic islands may be initiated by plasma compression at 
the shock front. Some synchrotron particles may be generated in our multiple magnetic 
reconnection process. Additionally, a phenomenon similar to striped wind may occur in the 
termination shock in our heliosphere H] • 
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